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We present a two-flavour linear sigma model with global chiral symmetry and (axial-)vector 
mesons as well as an additional glueball degree of freedom. We study the structure of the well- 
established scalar resonances /o(1370) and /o(1500): by a fit to experimentally known decay widths 
we find that /o(1370) is predominantly a qq state and /o(1500) is predominantly a glueball state. 
The overall phenomenology of these two resonances can be well described. Other assignments for 
our mixed quarkonium-glueball states are also tested, but turn out to be in worse agreement with 
the phenomenology. As a by-product of our analysis, the gluon condensate is determined. 
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I. INTRODUCTION 



Glueballs, the bound states of the gauge bosons of QCD, the gluons, are naturally expected in QCD due to the 
nonabelian nature of the theory: the gluons interact strongly with themselves and thus they can bind and form 
colorless states, analogously to what occurs in the quark sector. The existence of glueballs has been studied in the 
framework of the effective bag model for QCD already four decades ago [l| and it has been further investigated in a 
variety of approaches 0]. Numerical calculations of the Yang-Mills sector of QCD also find a full glueball spectrum 
D ■ in which the scalar glueball is the lightest state Q. 

Glueballs can mix with quarkonium {qq) states with the same quantum numbers. This makes the experimental 
search for glueballs more complicated, because physical resonances emerge as mixed states. The scalar sector J^'-^ = 
y—^ ' Q^'^ has been investigated in many works in the past. The resonance /o(1500) is relatively narrow when compared 
J> to other scalar-isoscalar states: for this reason it has been considered as a convincing candidate for a glueball state. 

Mixing scenarios in which two quark-antiquark isoscalar states fin and ss and one scalarglueball gg mix and generate 
^If!^ . the physical resonances /o(1370), /o(1500), and /o(1710) have been discussed in Refs. [JQ. 

' The aim of this work is to extend the linear chiral model of Refs. [1,0] by including the dilaton/glueball field. The 
• ] first attempt to incorporate a glueball into a linear sigma model was performed long ago in Ref. Q . The novel features 
' of the present study are the following: (i) The glueball is introduced as a dilaton field within a theoretical framework 
' where not only scalar and pseudoscalar mesons, but also vector and axial- vector mesons are present from the very 
I I beginning. This fact allows also for a calculation of decays into vector mesons. The model is explicitly evaluated 
• • ' for the case of Nf = 2, for which only one scalar-isoscalar quarkonium state exists: tr = fm which mixes with the 
, glueball. The two emerging mixed states are assigned to the resonances /o(1370) which is, in accordance with Ref. 
k> ' Q, predominantly a qq state, and with /o(1500) which is predominantly a glueball state, (ii) We consequently test 
; I , — to our knowledge for the first time — this mixing scenario above 1 GeV in the framework of a chiral model. 

' The model under consideration is built in accordance with the symmetries of the QCD Lagrangian. It possesses the 
known degrees of freedom of low-energy QCD [(pseudo)scalar and (axial-) vector mesons] as well as the same global 
chiral invariance. Another feature of the QCD Lagrangian is scale (or dilatation) invariance — \~^x^ (where is 
a Minkowski-space coordinate and A the scale parameter of the conformal group). It is realized at the classical level 
but broken at the quantum level due to the loop corrections in the Yang-Mills sector (scale anomaly). In this work 
the breaking of scale invariance is implemented at tree-level by means of a dilaton field (representing a glueball) with 
the usual logarithmic dilaton potential Q. However, all the other interaction terms (with the exception of the chiral 
anomaly) are dilatation-invariant in the chiral limit. 

Having constructed the Lagrangian of the effective model, we calculate the masses of the pure qq and glueball states 
in the J^^ = 0++ channel, study their mixing and calculate the decay widths of the mixed states. Although we work 
with Nf — 2, the use of flavor symmetry enables us to calculate the decay widths of the scalar resonances into kaons 
and into both the r\ and 77' mesons which contain the s-quark in their flavor wave functions. After the study of the 
already mentioned scenario where /o(1370) and /o(1500) are predominantly quarkonium and glueball, respectively, 
we also test the alternative scenario in which the resonance /o(1710) is predominantly glueball and scenarios in which 
/o(600) is predominantly quarkonium. They, however, lead to inconsistencies when compared to the present data and 
are therefore regarded as less favorable. 
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This paper is organized as follows. In Sec. II the Lagrangian of the linear sigma model with (axial-)vector and 
glueball degrees of freedom is constructed. In Sec. Ill we discuss the results for the masses of the quarkonium-glueball 
mixed states and their decay widths. In Sec. IV we present our conclusions. 



II. THE MODEL 



The Yang-Mills (YM) sector of QCD (QCD without quarks) is classically invariant under dilatations. This symmetry 
is, however, broken at the quantum level. The divergence of the corresponding current is the trace of the energy- 



momentum tensor TyJ^ of the YM Lagrangian 



YM) 



4g 



(1) 



where GJ^^ is the field-strength tensor of the gluon fields, g 



g{fj.) is the renornialized coupling constant at the scale 
fi, and the /3-function is given by f3{g) — dg/d\nfi. At the one-loop level /3{g) — —bg^ with b — 11Nc/{4:8tt'^). This 
implies g^{fJ,) = [2b\n{fj./AYM)]~ , where Aym — 200 MeV is the Yang-Mills scale. A finite energy scale thus emerges 
in a theory which is classically invariant under dilatation (dimensional transmutation). The expectation value of the 
trace anomaly does not vanish and represents the so-called gluon condensate: 



llNr 



where 



C"^ ~ (300 - 600 MeV)^ 



(2) 
(3) 

and lattice simu- 



The numerical values have been obtained through QCD sum rules (lower range of the interval) 
lations (higher range of the interval) (10| . 

At the composite level one can build an effective theory of the YM sector of QCD by introducing a scalar dilaton 
field G which describes the trace anomaly. The dilaton Lagrangian reads [8] 



1 TOq 

4 A2" 



G^n 
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(4) 



The minimum Gq of the dilaton potential is realized for Gq = A. Upon shifting G — >■ Go -|- G, a particle with mass 
ma emerges, which is interpreted as the scalar glueball. The numerical value has been evaluated in Lattice QCD and 
reads mc ^ 1-5 GeV (^j. The logarithmic term of the potential explicitly breaks the invariance under a dilatation 
transformation. The divergence of the corresponding current reads df_iJ^.^i — Tj^i ^ = —jItiqA^. This can be compared 

with the analogous quantity in Eq. ^ which implies A = •\/TTG^/(2mG). 

QCD with quarks is also classically invariant under dilatation transformations in the limit of zero quark masses 
(chiral limit). The scale of all hadronic phenomena is given by the previously introduced energy scale Aym- This 
fact holds true also when the small but nonzero values of the quark masses is considered. In order to describe these 
properties in a hadronic model we now extend the linear sigma model with U{Nf)R x U{Nf)L of Refs. by 
including the dilaton. To this end, the following criteria are applied [3| : (i) With the exception of the chiral anomaly, 
the parameter A from Eq. ([4]), which comes from the Yang-Mills sector of the theory in accordance with QCD, is 
the only dimensionful parameter of the Lagrangian in the chiral limit, (ii) The Lagrangian is required to be finite 
for every finite value of the gluon condensate Gq. This, in turn, also assures that no singular terms arise in the limit 
Go — > 0. In accordance with the requirements (i) and (ii) only terms with dimension exactly equal to 4 are allowed in 
the chiral limit. 

The hadronic Lagrangian obeying these requirements reads 



■'dil 



Tr 



(I?^$)t(i^^$)_™2|^J <i,tci, _ A2($t$)2 



Ai(Tr [$^$]; 



:[det($t) + det($)] -f^ Tr [H ($^ + $)] - ^Tr [{L^^'^f + {R^"'f] 



G 



2 j + ^^"^'^ + ^Tr[$tci,]Tr[L^L'^ + 

/i2Tr[$^L^L''$ + + 2/i3Tr[$i?^$^L^] + ... , 



(5) 
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where <J> represents the Nf x Nf (pseudo)scalar multiplet, and and i?^ the left- and right-handed vector multiplets. 
The dots represent further terms which do not affect the processes studied in this work and can be omitted. 

In the particular case oiNf = 2 studied in this manuscript one has $ = ((j + i77jv) + (ao + «7r) • t (our eta meson ryjv 



contains only non-strange degrees of freedom), L^^ — {u)^^ + /f ) i° + (/?" 



■ t and i?^ = (w^ - /C) t° + a'^) ■ t 



t are the generators of C/(2). Moreover, D^^ = ^/^^ ^i?/^), L^'' = L" - L^' , Rf^" = - d" . 

The explicit breaking of the global chiral symmetry is described by the term Tr[iJ($ + <i?^)] = ha [h — const. 
(X TTig), which allows us to take into account the nonvanishing value of the quark mass. This term contains the 
dimensionful parameter h with [h] = [energy^] and also explicitly breaks the dilatation invariance, just as the quark 
masses do in the underlying QCD Lagrangian. Finally, the chiral anomaly is described by the term c (det $ + det $t) 
[l5| . For Nf = 2 the parameter c carries the dimension [energy^] and represents a further breaking of dilatation 
invariance. This term, however, is driven by instantons, which are also a property of the Yang-Mills sector of QCD. 

The identification of the fields of the models with the resonances listed in Ref. [l^ is straightforward in the 
pseudoscalar and (axial-)vector sectors: the fields 7? and Ty^f correspond to the pion and the SU{2) counterpart of the 
r] meson, respectively, rjN = iuu + dd) / \/2, with a mass of about 700 MeV. This value can be obtained by 'unmixing' 
the physical t] and rj' mesons, which also contain ss contributions. The fields and represent the a;(782) 
and p(770) vector mesons, respectively, while the fields /C and represent the /i(1285) and ai(1260) axial- vector 
mesons, respectively. (In principle, the physical w and /i states also contain ss contributions but their admixture is 
small.) As shown in Ref. @, the a field should be interpreted as a predominantly qq state because its decay width 
decreases as 1/A"c in the limit of large number of colors. The tr and G fields mix: the physical fields a' and G' are 
obtained through an SO{2) rotation, as we shall show in the following. Then the first and most natural assignment is 
{cr',G'} = {/o(1370), /o(1500)}, see Sec. IIIA. Note that the oq state is assigned to the physical ao(1450) resonance 
in accordance with results of Ref. [Q]. Other assignments for {cr', G'} will be also tested in Sees. III.B and III.C and 
turn out to less favourable. 

In order to study the non-vanishing vacuum expectation values (vev's) of the two J^'-^ = 0++ scalar-isoscalar fields 
of the model a and G, we set all the other fields in Eq. dU to zero and obtain: 



G_ 

Go 



ha 



(6) 



Upon shiftingthe fields by their vev's, a — cr-f and G G+Go, we obtain the masses of the states a = (uu+dd) /^/2 
and G = 35 0, 



7^2 



Q2 ,ml^[l + i\u 



Gn 



A 



(7) 



Note that the pure glueball mass Mq depends also on the quark condensate 0, but correctly reduces to mc in the 
limit TOq = (decoupling of quarkonia and glueball). In the presence of quarkonia, mg 7^ 0, the vev Gq is given by 
the equation 



^^ = G;Jln 



Go 



A 



(8) 



The shift of the fields by their vev's introduces a bilinear mixing term ^ aG in the Lagrangian The physical 
fields a' and G' can be obtained through an 5*0(2) rotation. 



a' 
G' 



cos t/ sm t 
— sin 9 cos ( 



(9) 



with 



Ml, = Ml cos^ e + Ml SIT? + 2ml-^ sin(26') , 

Go 

Mg, = cos^ e + Ml sir? e-2ml-^ sin(26l) , 

Go 



(10) 
(11) 



where the mixing angle 9 reads 



■ arctan 



Go - M2 



(12) 
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The quantity ttiq can be calculated from the masses of the pion, rj^ and the bare a mass [see Ref. 



Z 



-Mi 



(13) 



If rriQ — c < 0, spontaneous breaking of chiral symmetry is realized. 



III. RESULTS AND DISCUSSION 



The Lagrangian ([5]) contains the following twelve free parameters: mg, Ai, A2, mi, gi, c, h, hi, /12, h^, togi 
A = vTTC^/(2tog). The processes that we will consider depend only on the combination hi + h2 + thus reducing 
the number of parameters to ten. We replace the set of ten parameters by the following equivalent set: 171^^, , nip, 
■niai, -^o-, WG, TO], C. The masses m^r (= 139.57 MeV) and mp {— 775.49 MeV) are fixed to their PDG values. 

As outlined in Refs. the mass of the t^n meson can be calculated using the mixing of strange and non-strange 

contributions in the physical fields -q and 77' (958): 

T] — rjN cos ip + Tjs sin ip, r]' — —rjN sin ip + r]s cos (p, (14) 

where ijs denotes a pure ss state and ip ~ —36° ^17]. In this way, we obtain the value m,,„ = 716 MeV. (Given the 
well-known uncertainty of the value of the angle (p, one could also consider other values, e.g., ip = —41.4°, as published 
by the KLOE Collaboration [31, which corresponds to = 755 MeV. Variations of the pseudoscalar mixing angle 
affect the results presented in this paper only slightly.) 

The value of m^i is fixed to 1050 MeV according to the study of Ref. JlS]- (We note that taking the present 
PDG estimate of 1230 MeV does not change our conclusions.) The chiral condensate is fixed as (p — Zf^^ and the 
renormalization constant Z is determined by the study of the process ai — >■ ttj: Z ~ 1.67 ± 0.2 @. 



A. Assigning a' and G' to /o(1370) and /o(1500) 

The cr' field denotes an isoscalar J^^ — 0++ state and its assignment to a physical state is a long-debated problem 
of low-energy QCD [1, 0, S E^- The two major candidates are the resonances /o(600) and the /o(1370) [16|. The 
study of Ref. @ has shown that /o(1370) is favoured to be a state which is predominantly qq. As stated above, 
the resonance /o(1500) is a convincing glueball candidate. For these reasons we first test the scenario in which 
{a', G'} = {/o(1370), /o(1500)}, which turns out to be phenomenologically successful, see below. 

We are left with the following four free parameters: C,Ma,raG,mi. They can be obtained by a fit to the five 
experimental quantities of Table |T| the masses of the resonances /o(1500) {Mqi = M^jj(i5oo) = 1505 MeV 16]) and 
/o(1370) (for which we use the mean value M^f = (1350 ± 150) MeV taking into account the PDG mass range 
between 1200 MeV and 1500 MeV [16]) and the three well-known decay widths of the well-measured resonance 
/o(1500): /o(1500) ^ tttt, /o(1500) ^ 7777 and /o(1500) -> KK. 



Quantity 


Our Value [MeV] 


Experiment [MeV] 




1191 ± 26 


1200-1500 


Mg' 


1505 ± 6 


1505 ± 6 


G' TTTT 


38 ±5 


38.04 ± 4.95 


G' TIT] 


5.3 ±1.3 


5.56 ± 1.34 


G' KK 


9.3 ±1.7 


9.37 ± 1.69 



TABLE I: Fit in the scenario {a , G'} = {/o(1370), /o(1500)}. Note that the /o(1370) mass ranges between 1200 MeV and 
1500 MeV [m S'Hd therefore, as an estimate, we are using the value ttIo-' = (1350 ± 150) MeV in the fit. 

Using the Lagrangian ([5]), these observables can be expressed as functions of the parameters listed above. Note 
that, although our framework is based on Nf = 2, we can calculate the amplitudes for the decays into mesons 
containing strange quarks by making use of the flavor symmetry SU (Nf — 3) [4] . It is then possible to calculate the 
following /o(1500) decay widths into pseudoscalar mesons containing s-quarks: /o(1500) KK, /o(1500) — > 7777 and 
/o(1500) ^ 7777'. 

The method yields xVd.o.f. = 0.29 (thus very smaU), C = (699 ± 40) MeV, = (1275 ± 30) MeV, ma = 
(1369 ± 26) MeV and mi = (809 ± 18) MeV. 
The consequences of this fit are the following: 
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(i) The quarkonium-glueball mixing angle reads 6 = (29.7 ± 3.6) °. This, in turn, imphes that the resonance /o(1500) 
consists to 76% of a glueball and to the remaining 24% of a quark-antiquark state. An inverted situation holds for 
/o(1370). 

(ii) Our fit allows us to determine the gluon condensate: C — (699 ± 40) MeV. This result implies that the upper 
value in Eq. ([2]) is favoured by our analysis. It is remarkable that insights into this basic quantity of QCD can be 
obtained from the PDG data on mesons. 

(iii) Further results for the /o(1500) meson are reported in the first two entries of Table HIl The decay into Air is 
calculated as a product of an intermediate pp decay. To this end the usual integration over the p spectral function is 
performed. Our result yields 30 MeV in the 47r decay channel and is about half of the experimental value rjg(i5oo)->.4,r 
= 54.0 ±7.1 MeV. However, it should be noted that an intermediate state consisting of two /o(600) mesons (which 
is also expected to contribute in this decay channel) is not included in the present model. The decay into the t]T]' 
channel is also evaluated; this channel is subtle because it is exactly on the threshold of the /o(1500) mass. Therefore, 
an integration over the spectral function of the decaying meson /o(1500) is necessary. The result is in a qualitative 
agreement with the experiment. 

(iv) The results for the /o(1370) meson are reported in the last four rows of Table ITTl They are in agreement with 
the experimental data regarding the full width: ^^(^ij^o) = (200-500) MeV [16||. Unfortunately, the experimental 
results in the different channels are not yet conclusive. Our theoretical results point towards a dominant direct tttt 
and a non-negligible 7777 contribution; these results correspond well to the experimental analysis of Ref. [20| where 
r/o(i370)^7r7r = 325 McV and ^ fo{i37o)^riri fo{i370)^TrTr = 0.19± 0.07 are obtained. We find that the four-pion decay 
of /o(1370) — > pp — > 47r is strongly suppressed (as was also found in Ref. Q). This is unlike Ref. [13], where a small 
but non-negligible value of about 50 MeV is found. However, it should be noted that due to interference effects our 
result for this decay channel varies strongly when the parameters are even slightly modified. 

(v) The mass of the p meson can be expressed as = ml + (jy^ {hi + h2 + /13) /2. In order that the contribution of 
the chiral condensate is not negative, the condition mi < mp should hold. In the framework of our fit this condition 
is fulfilled at the two-sigma level. This result points towards a dominant mi contribution to the p mass. This 
property, in turn, means that the p mass is predominantly generated from the gluon condensate and not from the 
chiral condensate. It is therefore expected that the p mass in the medium scales as the gluon condensate rather than 
as the chiral condensate. In view of the fact that mi is slightly larger than mp we have also repeated the fit by fixing 
mi = mp: the minimum has a x'^/d.o.f. ~ 1 and the results are very similar to the previous case. The corresponding 
discussion about the phenomenology is unchanged. 

(vi) As already stressed in Refs. the inclusion of (axial-) vector mesons plays a central role to obtain the present 
results. The artificial decoupling of (axial-)vector states would generate a by far too wide /o(1370) state. For this 
reason the glueball-quarkonium mixing scenario above 1 GeV has been previously studied only in phenomenological 
models with flavor symmetry IJ] but not in the context of chirally invariant models. 



Quantity 


Our Value [MeV] 


Experiment [MeV] 


G" pp — )• 4-K 


30 


54.0 ±7.1 


G' 7777' 


0.6 


2.1 ± 1.0 


a' TTTV 


284 ± 43 




o" — >■ 7777 


72 ±6 




a' -> KK 


4.6 ± 2.1 




a' ^ pp ^ Atv 


0.09 





TABLE II: Further results regarding the a' = /o(1370) and G' = /o(1500) decays. 



B. Assigning ct' and G" to /o(1370) and /o(1710) 

Although the resonance /o(1710) has also been regarded as a glueball candidate in a variety of works [2l|, its 
enhanced decay into kaons and its rather small decay width make it compatible with a dominant ss contribution in 
its wave function. Nonetheless, we have also tested the assumption that the pure quarkonium and glueball states mix 
to produce the resonances /o(1370) and /o(1710). 

The resonance /o(1710) is experimentally well known. Decays into tttt, KK, and 7777 have been seen, while no decay 
into 7777' and into An have been detected. Using the total decay width r^g(i7io) = (135 ± 8) MeV and the branching 
ratios reported in Ref. [l^ it is possible to deduce the decay widths into tttt, KK, and 7777, see Table III. 

A fit analogous to the one in Table I yields too large errors for the decay width a' = /o(1370) — >■ tttt. For this reason 
we repeat our fit by adding the following constraint: Ta'^Tnr = (250 ± 150) MeV. The large error assures that this 
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value is in agreement with experimental data on this decay width. The results of the fit are reported in Table III. 



Quantity 


Our Value [MeV] 


Experiment [MeV] 




1386 ± 134 


1350 ± 150 


Mg' 


1720 ± 6 


1720 ± 6 


G' -K-K 


29.7 ±6.5 


29.3 ±6.5 


G' nv 


6.9 ±5.8 


34.3 ± 17.6 


G' KK 


16± 14 


71.4 ±29.1 


a' TTTT 


379 ± 147 


250 ± 150 



TABLE III: Fit in the scenario {cr', G'} = {/o(1370), /o(1710)}. 



The fit yields xVd.o.f. = 1.72, C = (764 ± 256) MeV, M„ = (1516 ± 80) MeV, ma = (1531 ± 233) MeV and 
mi — (827 ± 36) MeV. The mixing angle between the pure quarkonium cr and the pure glueball G calculated from 
Eq. ([T2l) is 6* = 37.2° ± 21.4°. The is wors than in the previous case, but the overall agreement is acceptable. 
The mixing angle is large and could also overshoot the value of 45°, which would imply a somewhat unexpected and 
unnatural reversed ordering, in which /o(1370) is predominantly glueball and /o(1710) predominantly quarkonium. 

In Table IV we report the decay widths G" An, G' — ^ 7777', a' -> 7777 and a' — > KK , which can be calculated as a 
consequence of the fit of Table III. 



Decay Width 


Our Value [MeV] 


Experimental value [MeV] 


G' ^ 47r 


115 




G' ^ 


16 




a' — >■ 7777 


153 ± 79 




0' KK 


9 1+13.6 





TABLE IV: Further results from the fit with {a' , G'} = {/o(1370), /o(1710)}. 



A clear problem in the framework of this scenario emerges: the decay width G' = /o(1710) — 47r is large, while 
experimentally it has not been seen. Therefore, we conclude that this scenario is not favoured. Moreover, in this 
scenario the remaining resonance /o(1500) should then be interpreted as a predominantly ss state, contrary to what 
its experimentally dominant tttt decay pattern suggests. Consequently, /o(1710) is unlikely to be predominantly a 
glueball state; this is also in accordance with the results from the ZEUS Collaboration ^23i] . 

C. Scenarios with ct' = /o(600) 

The scenarios {cr',G"} = {/o(600), /o(1500)} and {cr',G"} = {/o(600), /o(1710)} have also been tested. In both 
cases the mixing angle turns out to be small (< 15°), thus the state /o(600) is predominantly quarkonium. Then, in 
these cases the analysis of Ref. applies: a simultaneous description of the tttt scattering lengths and the a — tttt 
decay width cannot be achieved. For these reasons the mixing scenarios with the resonance /o(600) as a quarkonium 
state are not favoured. 



IV. CONCLUSIONS AND OUTLOOK 

We have presented a globally chirally invariant linear sigma model with (axial-)vector mesons and a dilaton/glueball 
degree of freedom. We have studied the phenomenology of the scalar states for the case N j = 2: in the favoured 
scenario the resonance /o(1500) is predominantly a glueball with a subdominant qq component and, conversely, 
/o(1370) is predominantly a quark-antiquark [uu + dd)/^/2 state with a subdominant glueball contribution. It is 
interesting to observe that the success of the phenomenological description of these scalar resonances is due to the 
inclusion of the (axial-)vector mesons in the model. The gluon condensate is also an outcome of our study and turns 
out to be in agreement with lattice QCD results. Different scenarios in which /o(1710) is predominantly glueball 
and/or /o(600) is predominantly quarkonium do not seem to be in agreement with the present experimental data. 

Natural extensions of the model are the case Nf = 3 [l^] and the inclusion of a nonet of tetraquark states as 
additional low-lying scalar states. In this general scenario, a mixing of five scalar-isoscalar states takes place, which 
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allows to describe all relevant scalar-isoscalar resonances listed in the PDG below 1.8 GeV fl6'|. Applications of the 
model at nonzero temperature and density are also important because the presence of the dilaton field allows to study 
the restoration of both the dilatation and the chiral symmetry of QCD. 



Acknowledgments 

The authors acknowledge fruitful discussions with A. Heinz and E. Seel. The work of D.P. is supported by the 
Foundation " Polytechnical Society." 



[1] R. L. Jaffe and K. Johnson, Phys. Lett. B 60, 201 (1976); R. Konoplich and M. Shchepkin, Nuovo Cim. A 67, 211 (1982); 
M. Jezabek and J. Szwed, Acta Phys. Polon. B 14, 599 (1983); R. L. JafTe, K. Johnson and Z. Ryzak, Annals Phys. 
168, 344 (1986); M. Strohmeier-Presicek, T. Gutsche, R. Vinh Mau and A. Faessler, Phys. Rev. D 60, 054010 (1999) 
arXiv:hep-ph/9904461 . 

[2] F. E. Clo se, Rept. Prog. Phys. 51, 833 (1988); S. Godfrey and J. Napolitano, Rev. Mod. Phys. 71, 1411 (1999) 
[arXiv:hep -ph/9811410 ; C. Amsler and N. A. Tornqvist, Phys. Rept. 389, 61 (2004); E. Klempt and A. Zaitsev, Phys. 
Rept. 454, 1 (2007) arXiv:0708.4016 [hep-ph]]. 

[3] C. Morningstar and M. J. Peardon, AIP Conf. Proc. 688, 220 (2004) [arXiv:nucl-th/0309068]; M. Loan, X. Q. Luo and 
Z. H. Luo, Int. J. Mod. Phys. A 21, 2905 (2006) arXiv:hep-rat/0503038 ; R^. Gregory, A. C. Irving, C. C. McNeile, 
S. Miller and Z. Sroczynski, PoS LAT2005, 027 (2006) arXiv:hep-lat/0510066 ; Y. Chen et al, Phys. Rev. D 73, 014516 
(2006) iarXiv:h ep-lat/0510074 . 

[4] F. Giacosa, T. Gutsche, V. E. Lyubovitskij and A. Faessler, Phys. Rev. D 72, 094006 (2 005) |arX iv:hep-ph/05 09247| ; 

F. Giacosa, T. Gutsche, V. E. Lyubovitskij and A. Faessler Phys. Lett. B 622, 277 (2005) [arXiv:he p-ph/0504033 . 
[5] C. Amsler and F. E. Close, Phys. Rev. D 53, 295 (1996) arXiv:hep-ph/9507326 ; F. E. Close and A. Kirk, Eur. Phys. 

J. C 21, 531 (2001) arXiv:hep-ph/0103173 ; F. Giacosa, T. Gutsche and A. Faessler, Phys. Rev. C 71, 025202 (2005) 

arXiv:hep-ph/0408085 ; 

[6] D. Parganhja, F. Giacosa and D. H. Rischke, Phys. Rev. D 82, 054024 (2010) arXiv:1003.4934l [hep-ph]]. 

[7] S. Gallas, F. Giacosa and D. H. Rischke, Phys. Rev. D 82, 014004 (2010) arXiv:0907.5084 [hep-ph]]; D. Parganlija, 

F. Giacosa and D. H. Rischke, Acta Phys. Polon. Supp. 3 (2010) 963 arXiv: 1004.4817 [hep-ph]]. 

[8] C. Rosenzweig, A. Salomone and J. Schechter, Phys. Rev. D 24, 2545 (1981); A. Salomone, J. Schechter and T. Tudron, 
Phys. Rev. D 23, 1143 (1981); C. Rosenzweig, A. Salomone and J. Schechter, Nucl. Phys. B 206, 12 (1982) [Erratum-ibid. 
B 207, 546 (1982)]; A. A. Migdal and M. A. Shifman, Phys. Lett. B 114, 445 (1982); H. Gomm and J. Schechter, Phys. 
Lett. B 158, 449 (1985); R. Gomm, P. Jain, R. Johnson and J. Schechter, Phys. Rev. D 33, 801 (1986). 

[9] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979); L. J. Reinders, H. R. Rubinstein and 
S. Yazaki, Nucl. Phys. B 186, 109 (1981); J. Marrow, J. Parker and G. Shaw, Z. Phys. C 37, 103 (1987); B. V. Geshkenbein, 
Sov. J. Nucl. Phys. 51, 719 (1990) [Yad. Fiz. 51, 1121 (1990)]; D. J. Broadhurst, P. A. Baikov, V. A. Ilyin, J. Fleischer, 
O. V. Tarasov and V. A. Smirnov, Phys. Lett. B 329, 103 (1994) arXiv:hep-ph/9403274 ; F. J. Yndurain, Phys. Rept. 320, 
287 (1999) arXiv:hep-ph/9903457 ; B. L. loffe and K. N. Zyablyuk, Eur. Phys. J. C 27, 229 (2003) ar Xiv:hep-ph/0207183) ; 
K. Zyablyuk, JHEP 0301, 081 (2003) arXiv;hep-ph/0210103 ; A. Samsonov, arXiv;hep-ph/0407199 
[10] J. Kripfganz, Phys. Lett. B 101, 169 (1981); A. Di Giacomo and G. C. Rossi, Phys. Lett. B 100, 481 (1981); A. Di 
Giacomo and G. Paffuti, Phys. Lett. B 108, 327 (1982); E. M. Ilgenfrhz and M. MuUer-Preussker, Phys. Lett. B 119, 395 
(1982); S. s. Xue, Phys. Lett. B 191, 147 (1987); M. Campostrini, A. Di Giacomo and Y. Gunduc, P hys. Lett. B 225, 393 
(1989); A. Di Giacomo, H. Panagopoulos and E. Vicari, Nucl. Phys. B 338, 294 (1990); X. D. Ji, |arXiv:hep-ph /9506413| 

G. Boyd and D. E. Miller, arXiv:hep-ph/9608482 

[11] S. Gasiorowicz and D. A. Geffen, Rev. Mod. Phys. 41, 531 (1969); J. Boguta, Phys. Lett. B 120, 34 (1983); P. Ko and 

S. Rudaz, Phys. Rev. D 50, 6877 (1994). 

[12] D. Parganlija, F. Giacosa and D. H. Rischke, PoS CONFINEMENT8, 070 (2008) arXiv:0812 . 21831 [hep-ph]]. 
[13] M. Urban, M. Buballa and J. Wambach, Nucl. Phys. A 697, 338 (2002) .arXiv:hep-ph/0102260| . 
[14] F. Giacosa, Phys. Rev. D 80, 074028 (2009) arXiv:0903.4481 [hep-ph]]. 
[15] G. 't Hooft, Phys. Rept. 142, 357 (1986). 

[16] K. Nakamur a etal. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

[17] F. Giacosa. larXiVr07 12.0186 [hep-ph]. 

[18] F. Ambrosino et ai, JHEP 0907, 105 (2009) 'arXiv:0906.3819 [hep-ph]]. 

[19] R. L. Jaffe, Phys. Rev. D 15, 267 (1977); R. L. Jaffe, Phys. Rev. D 15, 281 (1977); A. H. Fariborz, R. Jora and J. Schechter, 
Phys. Rev. D 72, 034001 (2005); A. H. Fariborz, Int. J. Mod. Phys. A 19, 2095 (2004); M. Napsuciale and S. Rodriguez, 
Phys. Rev. D 70, 094043 (2004); F. Giacosa, Phys. Rev. D 75, 054007 (2007). 

[20] D. V. Bugg, Eur. Phys. J. C 52, 55 (2007) arXiv:0706.1341 [hep-ex]]. 

[21] W. J. Lee and D. Weingarten, Phys. Rev. D 61, 014015 (2000) ;arXiv:hep-lat/9910008| ; H. Y. Cheng, C. K. Chua and 
K. F. Liu, Phys. Rev. D 74, 094005 (2006) arXiv:hep-ph /0607206| " 



8 



[22] D. Parganlija, F. Giacosa and D. H. Rischke, arXiv:0911.3996 [nucl-th]. 

[23] S. Chekanov et al. [ZEUS Collaboration] Phys. Rev. Lett . 101, 112003 (2008) [arXiv: 0806. 08071 [hep-ex]] : U. Karshon [HI 
Collaboration and ZEUS Collaboration], larXiv:0907.4574l [hep-ex]. 

[24] D. Parganlija, F. Giacosa, D. H. Rischke, P. Kovacs and G. Wolf, Int. J. Mod. Phys. A 26, 607 (2011) >rXiv:1009.2250 [hep- 
ph]]; D. Parganlija, F. Giacosa, P. Kovacs and G. Wolf, larXiv : 1 1 iT6 1 4 [hep-phj; D. Parganlija, F. Giacosa, D. H. Rischke, 
P. Kovacs and G. Wolf, in preparation. 



